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Short  Title  of  Work 


Pre-Breakdown  Laser  Target  Vaporization  and 
Enhanced  Thermal  Coupling 


ABSTRACT 


\ aporization  of  metallic  targets  under  the  influence  of  pulsed 
C°2  laser  radiation  is  observed  prior  to  breakdown  plasma  formation 
using  the  techniques  of  atomic  absorption  spectroscopy  and  time 
resolved  emission  spectroscopy.  Aluminum  targets  vaporize  in  times 
short  compared  to  predictions  based  on  bulk  thermal  properties,  a 
behavior  which  may  be  attributable  to  microsurface  imperfections. 
Plasma  formation  times  are  consistent  with  models  that  assume  the 
presence  of  metallic  vapor. 

Aluminum  targets  that  experienced  plasma  breakdown  exhibited 
cratering  vs  little  visible  surface  damage  to  targets  irradiated  at 
fluxes  just  below  the  breakdown  threshold.  Thermal  energy  deposition 
per  unit  incident  energy  was  increased  by  a factor  of  three  when  the 
breakdown  threshold  was  exceeded. 
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I.  INTRODUCTION 
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Laser  supported  absorption  waves  are  a well  documented 
phenomenon  in  which  a laser  beam's  energy  is  partially  or  fully  absorbed 
in  a front  of  ionized  gas  which  is  maintained  by  the  absorption  process. 
These  waves  fall  into  two  broad  categories.  Laser  supported  detonation 
(LSD)  waves  are  characterized  by  absorption  lengths  less  than  the  laser 
beam  diameter,  pressures  considerably  greater  than  ambient  pressure, 
and  supersonic  velocities.  Laser  supported  combustion  waves  (LSC) 
attenuate  the  beam's  energy  over  a distance  on  the  order  of  or  greater 
than  the  beam  diameter,  the  pressure  is  on  the  order  of  ambient  and 

the  velocities  are  less  than  sonic. 

It  is  also  well  documented  both  theoretically  and  experimentally 
that  breakdown,  i.e.  , the  creation  of  substantial  ionization,  in  clean 
air  at  STP  requires  laser  fluxes  of  greater  than  2 x 10^  watts /cm  for 
10.  6 jjl  wavelength  radiation.  Under  such  intensities  a detonation  wave 
results.  However,  it  has  been  observed  that  breakdown  and  subsequent 
wave  propagation  can  occur  at  substantially  lower  fluxes  if  the  beam 
impinges  on  a target  consisting  of  any  one  of  a broad  range  of  materials. 

The  ignition  of  absorption  waves  can  have  a number  of  important 
implications  with  regard  to  laser  applications.  It  is  clear  that  a highly 
ionized  absorption  wave  that  propagates  a sufficient  distance  away  from 
a target  will  shield  the  target  from  further  absorption  of  radiation.  If 
the  wave  is  in  detonation  regime,  the  high  pressure  volume  recedes 
from  the  target  leading  to  a decrease  in  the  impulsive  coupling. 
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On  the  other  hand,  there  is  experimental  evidence  indicating  that 
a target  with  small  absorptivity  to  the  incoming  radiation  can  experience 
enhanced  thermal  coupling  if  an  absorption  wave  is  ignited.  Possible 
explanations  for  this  effect  are  re-radiation  by  the  highly  absorbing 
plasma  back  to  the  target  and/or  convective,  conductive  transfer  of  heat 
to  the  target.  There  are  also  theoretical  predictions  that  a detonation 
wave  which  is  allowed  to  propagate  no  further  than  one  laser  beam 
diameter  away  from  the  target  can  efficiently  couple  impulses  to  the 
target.  ^ 

Since  target  effects  can  be  dramatically  dependent  on  the  presence 
or  non-presence  of  a laser  supported  absorption  wave,  it  is  important 
to  know  under  what  conditions  - target  material,  surface  temperature 
and  roughness;  ambient  gas  pressure  and  flow;  laser  flux,  pulse  time 
and  spot  size  - these  waves  are  ignited.  A fair  body  of  empirical  data 
already  exists.  However,  lacking  a suitably  tested  model  for  the 
mechanism  responsible  for  the  ignition  of  these  waves,  the  present  data 
has  not  been  correlated  ror  is  it  possible  to  extrapolate  to  conditions 
beyond  those  tested.  This  latter  point  is  particularly  important  if  one 
is  to  project  possible  countermeasures  either  offensive  or  defensive. 

It  is  the  purpose  of  the  present  investigation  to  determine  the 
mechanism  of  wave  ignition  unde  a limited  set  of  conditions  which  are 
of  interest  to  fviture  applications  of  pulsed  laser  devices.  Specifically, 
this  investigation  will  deal  with  fluxes  in  the  range  of  10^  to  a few  times 
10  watts/cm  at  the  largest  spot  sizes  consistent  with  the  10^  watt  10.  6 fj. 
laser  output.  Most  of  the  work  will  deal  with  metallic  targets. 
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II.  BREAKDOWN  INITIATION  STUDIES 


Three  principal  mechanisms  have  been  proposed  to  explain  surface 

8 2 

plasma  breakdown.  At  high  incident  flux  values  (>  10  watts/cm  ) the  elec- 
trie  field  associated  with  the  wave  would  be  sufficient  to  induce  field  emis- 
sion from  needle  shaped  points  on  the  material  surface.  Electrons  would 
gain  sufficient  energy  from  the  field  within  their  mean  free  path  to  ionize 
an  air  molecule.  A cascade  process  follows  leading  to  full  ionization.  A 
s milar  model,  at  somewhat  lower  fluxes,  derives  the  initial  electrons  from 
thermionic  emission.  Howe  er,  it  seems  likely  that  points  hot  enough  to 
emit  thermionically  would  soon  vaporize.  Not  only  would  this  remove  the 
emitting  points,  but  the  resulting  va~'  'r  would  dominate  the  subsequent 
physics . 

b 8 2 ^ / 

At  fluxes  of  10  to  10  watts/cm  most  theoretical  investigations  ’ 
have  assumed  that  the  target  first  vaporizes  and  that  the  initial  electron 
concentration,  derivable  from  the  Saha  equation  evaluated  at  the  vaporiza- 
tion temperature,  then  increases  due  to  a cascade  ionization  process  in  tk** 
ta rget  vapor. 

Although  these  models  are  reasonably  successful  in  predicting  cas- 
cade ionization  times  in  the  vapor  which  are  consistent  with  observed  break- 
down times,  they  all  must  call  for  anomalous  target  absorptivities  or 
thermal  conductivities  in  order  to  generate  target  vapor  on  time  scales  con- 
sistent with  observed  breakdown  times.  For  instance,  an  Al  surface  is  ex- 

z.  2 

pected  to  vaporize  in  27  gsec  with  an  absorbed  flux  of  10  watts  per  cm  . 


However,  the  initial  absorptivity  is  » 3%  which  implies  27  Msec  for  incident 

fluxes  of  3 x 10  watts  cm2.  Breakdown  is  obse  rved  on  a time  scale  of 
/i  secs  at  this  flux. 

One  explanation  for  rapid  vaporization  would  be  heating  of  small 
flakes  or  points  on  the  surface  which  are  in  poor  thermal  contact  with  the 
bulk  material.  The  fact,  observed  in  our  studies  as  well  as  others,  that  a 
surface  that  does  not  break  down  at  a given  flux  can  be  induced  to  do  so  by 
roughening  it  with  emery  cloth  gives  substance  to  this  view. 

One  of  the  principal  objectives  of  this  investigation  was  to  determine 
if  target  vapor  actually  evolves  from  the  target  previous  to  breakdown,  re- 
gardless of  the  surface  mic rostructu ral  details  that  lead  to  the  vaporiza- 
tion. Detection  of  the  vapor  would  greatly  enhance  the  credence  one  could 
place  on  the  theoretical  models  that  assume  vaporization. 

Since  the  required  vapor  for  breakdown  can  be  quite  minute  and 
since  the  vapor  need  not  be  strongly  optically  emitting  at  its  vaporization 
temperature,  a particularly  sensitive  technique  is  needed  to  detect  its  pres- 
ence. To  this  end  we  have  employed  a modification  of  the  method  of  Res- 
onance Absorption  Spectroscopy7  which  has  long  been  used  by  chemists  for 
detecting  trace  elements. 

A.  ATOMIC  ABSORPTION  SPECTROSCOPY 

1 . Theoretical  Background 

The  absorption  coefficient  for  photons  with  a frequency  correspond- 
ing to  the  energy  difference  between  two  atomic  states  is  given  by  the  ex- 
pression 

tt  c r f 

a = N - - e 

A v 
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where  a is  the  absorption  coefficient  in  inverse  cm,  N is  the  number  den- 
sity of  atoms  in  the  lower  state  of  the  transition,  rf  is  the  classical  elec- 
tron radius,  c is  the  speed  of  light  and  f is  the  oscillator  strength  for  the 
transition  involved.  The  width  of  the  transition.  Av  in  frequency  units  is 
governed  by  a number  of  considerations.  At  low  densities  in  a cold  gas, 

A v is  the  natural  line  width  which  is  given  by  the  inverse  of  the  lifettme  of 
the  transition.  In  the  limit  of  high  density,  resonance  broadening,  in  which 
excitation  from  an  excited  atom  is  transferred  to  an  unexcited  atom  of  the 
same  species,  dominates  and  the  frequency  width  is  given  by 


v . i v!  Nf  (iiV/2 
2 v w 


where  g.  and  gu  are  the  multiplicities  of  the  lower  and  upper  states  respec- 
l u 

tively.  The  absorption  coefficient  in  this  case  is 

2 pu\l/2i 

" = 2 ' w x 

that  is,  the  radiation  is  absorbed  over  a distance  on  the  order  of  a wave- 
length, A.  independent  of  the  density  or  the  f number.  This  condition  of 
high  density  could  be  expected  tmmediately  above  the  surface  of  a strongly 

vaporizing  target. 

The  most  generally  expected  broadening  mechanism  over  the  range 
of  experimental  conditions  would  be  due  to  the  thermal  motions  of  the  hot 
vapor  giving  a Doppler  width  of 

Av  = y2  7rkTa/Ma"  v/c 
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where  T and  M are  respectively  the  temperature  and  atomic  mass  of 
a a 

species  a.  This  mechanism  is  dominant  for  densities  satisfying 

1/2 


N « 1 


where  0 is  the  vapor  temperature  in  eV,  A&  is  the  atomic  mass  number 
3. 

and  X is  the  wavelength  of  the  transition  in  fi . 

4 

As  a numerical  example,  we  consider  aluminum  vapor  at  the  vapor- 
ization temperature  of  2750°K  and  a ground  state  population  of  N = 3.  6 x 
10U  /cm3.  Virtually  all  the  electrons  will  be  in  the  ground  state  so  N also 
represents  the  total  atomic  density.  Choosing  for  example  the  aluminum 
resonance  transition  at  X = 3092  R with  an  f number  of  0.79,  we  find  that 
Doppler  broadening  is  dominant  and  the  absorption  coefficient  is  l/cm. 

This  number  illustrates  the  extreme  sensitivity  of  the  method  to  trace 
quantities  of  target  vapor. 

2.  Experimental  Method 

The  basic  experimental  plan  has  two  configurations.  An  appropriate 
source  of  resonance  radiation  light  for  the  desired  target  vapor  transition 
is  either  directed  across  the  target,  as  shown  in  Fig.  1.  or  else  is  aimed 
to  shine  on  the  same  target  spot  as  the  impinging  CC>2  laser  radiation,  as 
shown  in  Fig.  2.  The  light  passing  by  the  target  or  reflecting  from  it  is 
collected  with  a lens  and  directed  to  a Jarrell  Ash  0.  25  meter  mono- 
chrometer which  serves  to  filter  out  undesired  lines  from  the  resonance 
radiation  source  and  to  partially  block  light  due  to  a developing  target 
breakdown.  The  transmitted  radiation  is  detected  with  a 1P28  photomulti- 
plier. 
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ATOMIC  ABSORPTION  SPECTROSCOPY 
(TRANSVERSE  LIGHT) 


Fig.  1 Schematic  of  Experimental  Apparatus 


ATOMIC  ABSORPTION  SPECTROSCOPY 
(REFLECTED  LIGHT) 


In  the  first  configuration,  the  spatial  resolutio.  (measured  experi- 
mentally by  observing  the  intensity  of  transmitted  light  as  a function  of  dis- 
tance as  a sharp  edged  blade  is  moved  across  the  focal  spot  near  the  target) 
is  - 1 mm  as  shown  in  Fig,  3.  This  leads  to  some  ambiguity  in  the  density 
of  the  target  vapor  since  partial  transmittal  of  light  could  be  due  to  very 
dense  completely  absorbing  vapor  that  has  moved  only  part  way  across  the 
resolution  element.  On  the  other  hand,  this  configuration  has  the  advantage 
that  the  probe  beam  can  be  positioned  at  various  distances  from  the  target 
allowing  the  velocity  of  the  vapor  to  be  obtained  by  noting  time  of  onset  of 
probe  beam  absorption.  This  information  can  in  turn  be  used  to  offset  the 
aforementioned  ambiguity. 

If  a target  were  to  vaporize  uniformly,  all  the  probe  light  in  the  sec- 
ond configuration  would  intercept  the  vapor  as  soon  as  it  left  the  target. 
However,  this  method  suffers  several  disadvantages.  Firstly,  changes  in 
the  target  reflectivity  or  scattering  coefficient,  certainly  expected  in  the 
case  of  prevaporization  melt,  would  lead  to  changes  in  the  received  signal. 
Secondly,  the  vapor  has  a velocity  component  in  the  direction  of  the  probe 
beam  which  shifts  the  center  of  the  absorption  line  away  from  the  center  of 
the  probe  line  radiation.  For  these  reasons,  only  limited  data  were  taken 
by  this  method  and  they  will  not  be  discussed  here. 

The  resonant  radiation  employed  in  these  experiments  was  obtained 
from  standard  hollow  cathode  lamps  manufactured  by  Perkin  Elmer  Corp. 
for  use  in  flame  spectroscopy  apparatus.  Different  lamps  were  obtained 
with  cathodes  constructed  of  the  same  material  as  the  test  target  material 
under  investigation.  These  materials  included  Zn,  Al  and  Ti.  The  lamps 
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distance  from  target 


Transmitted  intensity  as  a Faction  oi  Biade  Position  in  Atomic 
Absorption  Spectroscopy  Probe  Beam 


contain  a noble  gas,  usually  Ne  or  Ar,  and  the  emitted  radiatic  n consists 
of  lines  from  these  gases  as  well  as  lines  from  low  density  vaporized 
cathode  material. 

The  standard  operating  current  for  these  lamps  is  on  the  order  of 
10  mA  which  is  more  than  adequate  for  their  usual  application.  However, 
it  was  found  that  the  intensity  of  the  emitted  radiation  was  insufficient  for 
giving  acceptable  signal-to-noise  ratios  when  time  resolution  on  the  order 
of  1 gisec  is  demanded  as  is  the  case  in  this  investigation.  To  correct  this 
situation,  the  lamps  were  maintained  at  their  normal  10  mA  DC  level  and 
then  pulsed  to  a 2 amp  level  for  a 200  p.sec  interval  which  included  the  time 
of  the  laser  target  interaction. 

The  "Breakdown  Laser",  an  E-beam  sustainer  CC^  laser  device, 
was  used  to  irradiate  samples  of  Zn,  Ti  and  Al.  The  pulse  length  was 
100  fj.  sec  with  energies  ranging  up  to  100  J and  the  diameter  of  the  first 
Airy  ring  was  4 mm.  The  laser  energy  was  adjusted  for  a given  sample  so 
as  to  be  just  below  or  just  above  the  breakdown  level.  Fluxes  exceeding  the 
breakdown  threshold  by  even  20%  would  cause  breakdown  in  times  too  short 
(less  than  a few  gsec)  to  resolve  the  onset  of  target  vaporization  before  the 
atomic  absorption  apparatus  was  overwhelmed  with  breakdown  light. 

3.  Experimental  Results 

Figure  4 shows  typical  results  obtained  with  a Ti  target  at  a flux  just 
below  breakdown.  The  radiation  from  a titanium  hollow  cathode  lamp  passes 
within  1 mm  of  the  target  surface.  The  presence  of  titanium  vapor  is  sig- 
naled by  the  onset  of  absorption  of  the  3643  J?  radiation  within  ~ 8 jusec  of 
the  laser  turn  on.  As  a control,  radiation  from  a zinc  hollow  cathode  lamp 
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03  20/x  SEC / DIV ► 

Fig.  4 (a)  Transmission  of  3643  R Ti  Light  above  a Ti  Target 

(b)  Laser  Output 


was  also  passed  above  the  target  on  subsequent  shots  (see  Fig.  5).  The 
fact  that  the  nearby  wavelength,  3076  R,  radiation  is  not  affected  shows 
that  the  decreases  in  signal  observed  with  the  Ti  lamp  are  not  due  to  dis- 
toi  tiono  of  the  optical  medium,  but  rather  to  resonance  absorption  by  the 
Tl  vaPor-  The  experiments  were  also  repeated  with  an  A1  lamp  at  3093  R 
(Fig.  6).  Here  even  stronger  absorption  is  observed  and  at  an  earlier 
time  t ~ 4 M sec,  after  laser  turn  on.  The  alloy  employed  was  Ti  6,  4 
consisting  of  90%  Ti,  6%  A1  and  4%  vanadium.  It  appears  that  the  aluminum 

constituent  is  distilling  earlier  from  the  sample  due  to  its  lower  boiling 
point. 

Experimental  results  for  A1  targets  are  presented  in  Fig.  7.  The 
outputs  from  the  geld  doped  germanium  laser  power  monitor,  the  lower 
trace  in  each  oscillogram,  are  unfiltered  and  show  peaks  of  instantaneous 
intensity  on  the  order  of  twice  the  average  on-axis  intensity  of  ^ 7 x 106 
watts/cm  . The  spikiness  of  the  laser  output  is  due  to  parasitic  oscilla- 
tions induced  by  feedback  from  the  target.  The  top  traces  in  each  oscillo- 
gram show  the  transmission  f the  probe  light  and  the  zero  transmission 
reference  level.  In  oscillogram  (a)  absorption  of  resonant  A1  3093  R light 
appears  within  about  3 Msec  of  laser  turn  on.  The  absorption  increases  to 
approximately  60%  at  a time  30  gsec  after  laser  turn  on.  The  trace  then 

disappears  due  to  the  large  output  of  light  associated  with  the  plasma  break- 
down that  occurs  at  that  time. 

After  breakdown,  the  target  spot  appeared  shiny  against  the  as- 
received  rolled  surface.  A following  section  will  discuss  the  morphology 
of  this  melt  zone  in  some  detail.  If  the  same  spot  was  hit  again  with  the 
laser,  the  results  shown  in  oscillogram  (b)  were  obtained.  There  was  no 
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evidence  of  target  vaporization  and  no  plasma  breakdown  ensued.  Subse- 
quent shots  would  not  ignite  a breakdown  even  in  the  presence  of  intention- 
ally introduced  surface  contamination  in  the  form  of  fingerprints.  However, 
a slight  roughening  of  the  surface  with  carborundum  paper  would  ensure  a 
breakdown  on  the  following  shot. 

With  laser  power  slightly  above  the  threshold  for  breakdown,  the 
time  to  breakdown  varied  randomly  from  a microsecond  to  tens  of  ^ sec . 
Absorption  of  A1  resonant  radiation  was  seen  in  all  cases  for  which  there 
was  sufficient  time  for  the  vapor  to  block  a significant  portion  of  the  probe 
beam  and  which  were  within  the  time  response  of  the  instrumentation.  For 
those  shots  in  which  breakdown  occurred  within  2 fj.sec  or  less  after  laser 
turn  on,  such  as  shown  in  oscillogram  (c),  no  direct  evidence  of  vaporiza- 
tion was  obtained.  There  is,  however,  no  reason  to  assume  that  vapor  was 
not  present  and  that  a new  breakdown  mechanism  was  operative  in  those 
cases. 

Oscillogram  (d)  shows  the  transmission  above  the  Al  target  of 
3076  R light  from  a Zn  lamp.  Here  there  is  no  absorption  before  the  trace 
goes  off  scale  at  the  time  of  breakdown,  once  again  proving  that  true  res- 
onant absorption  as  opposed  to  refraction  or  scattering  is  being  measured. 

One  case  of  particular  interest  was  Zn,  for  which  vapor  was  ob- 
served out  to  a distance  of  several  cm  in  front  of  the  target  for  times  up  to 
1/2  sec  after  breakdown.  Figure  8 shows  the  increase  of  the  transmitted 
Zn  lamp  signal  on  a slow  100  msec/div  time  scale.  The  irregular  nature  of 
the  increase  was  due  to  the  random  effects  of  the  air  currents  that  con- 
vected  the  vapor  laden  air  out  of  the  observation  region.  It  should  be  noted 
that  the  Zn  vapor  was  by  this  time  quite  cold  and  nonemitting.  Indeed, 
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'ig.  8 (a,  c)  Transmission  of  3076  $ Zn  Light  above  a Zn  Target 

(b)  Laser  Output 


interferograms  of  the  breakdown,  Fig.  9,  show  that  the  medium  has  re- 
amed to  a -elatively  undisturbed  condition  by  10  msec  after  the  laser 
pulse.  These  results  suggest  that  caution  should  be  used  in  interpreting 
multiple  pulse  experiments  in  which  the  air  in  front  of  the  target  is  not 
changed  between  shots.  Low  levels  of  vapor  ur detectable  by  other  means, 
could  linger  in  the  target  area  and  trigger  breakdowns  away  from  the  tar- 
get surface. 

B,  TIME  RESOLVED  EMISSION  SPECTROSCOPY 

Atomic  absorption  spectroscopy  provides  a sensitive  probe  for  the 
detection  of  most  metallic  vapors.  The  method  does  have  its  limitations. 
As  mentioned,  it  is  difficult  to  "see"  vapor  that  is  near  the  target  surface. 
Moreover,  the  number  of  sources  of  resonance  radiation  is  limited  and  it 
is  quite  possible  that  a surface  contaminant,  which  would  be  undetected  in 
a given  experiment,  could  be  evolving  from  the  target  and  contributing  to 
the  breakdown  process.  Emission  spectroscopy  covering  a broad  spectral 
range  can  be  useful  here,  given  sufficient  sensitivity  and  time  resolution. 
Unfortunately,  it  may  bi  shown  that  ordinary  photographic  film  does  not 
have  the  sensitivity  to  record  emission  from  spectral  lines  radiating  at  in- 
tensities corresponding  to  a black  body  at  A1  vaporization  temperature  if 
time  resolution  of  jj.secs  are  demanded  and  spectrometers  of  reasonable  f 
numbers  are  assumed. 

1 . Experimental  Method 

In  answer  to  these  shortcomings,  a spectroscopic  facility  was  as- 
sembled in  which  a gateable  Silicon  Intensifier  Tube  (S.  I.  T.  tube)  Optical 
Multichannel  Analyzer  (O.  M.  A.  ) manufactured  by  SSR  Instrument  Co.  was 
used  to  record  time  resolved  spectra.  In  essence,  a visible  spectrogram 
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Fig.  9 Breakdown  on  an  Unpolished  Zinc  Target 
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from  a grating  spectrograph  is  projected  on  the  face  of  the  S.  I T photo- 
cathode. Photoelectrons  from  the  cathode,  whose  passage  is  gateable  to 
times  as  short  as  100  nsec  in  synchronism  with  the  experiment,  are  ac- 
celerated to  a silicon  target  where  a charge  image  is  formed  at  a gain  of 
2000.  This  one-dimensional  charge  image  is  then  read  by  a cathode  ray 
beam  and  the  resulting  information  is  stored  digitally  in  a 500  element 
memory  bank.  The  data  may  be  read  out  at  any  time  to  plot  a complete 
spectrogram,  or  the  individual  spectral  lines  may  be  displayed  and  the 
area  under  them  integrated. 

A Jarrell  Ash  1 /4  meter  monochromator  with  a 300  line  per  mm 
grating  was  modified  for  mounting  the  S.  I.  T.  tube  at  the  slitless  exit 
plane.  This  provided  a 1500  ^ spectral  coverage  at  3 ^ per  channel  with 
the  central  wavelength  variable  over  the  tube  response  range  (S20  photo- 
cathode) by  a grating  rotation. 

Successful  application  of  the  device  demanded  elimination  of  sev- 
eral sources  of  interference.  Electrical  pickup  was  avoided  by  slaving 
the  laser  to  fire  at  a particular  time  with  respect  to  the  internal  clock 
cycle  of  the  O.  M.  A.  unit.  TheS.I.T.  tube  proved  to  be  quite  sensitive 
to  microphonics.  This  was  a severe  problem  since  shock  waves  result- 
ing from  a laser  produced  plasma  breakdown  disturbed  the  tube  during 
the  time  that  the  silicon  charge  storage  element  was  being  read  out  b/ 
the  electron  beam.  Gross  noise  in  the  spectra  was  finally  removed  by 
mounting  the  spectrometer  tube  assembly  in  a lead  box  lined  with  foam. 
The  lens  used  to  focus  the  target  on  the  entrance  slit  served  as  a window 
for  the  box.  These  precautions  provided  the  requisite  acoustical  damping. 
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A solute  intensity  calibration  of  the  unit  was  obtained  by  placing  a 
standard  tungsten  lamp  such  that  its  filament  occupied  the  position  in  the 
optical  train  that  would  normally  be  occupied  by  the  laser  target.  Calibra- 
tion was  performed  with  the  S.  I.T.  tube  gated  on  for  the  same  times  as 
were  used  for  recording  plasma  breakdown  spectra.  An  electro-mechanical 
shutter  of  20  ms  open  duration  was  placed  between  the  filament  and  entrance 
slit  to  eliminate  effects  of  imperfect  blockage  of  light  when  the  S.I.  T.  tube 
was  in  its  off  state. 

2.  Experimental  Results 

Speccral  emissions  from  carbon  targets  were  observed  at  an  angle 
of  45°  from  the  target  normal.  The  gate  pulse  to  the  S.  I.  T.  tube  was  set  at 
1 ju  sec  width  and  the  delay  between  the  initiation  of  the  laser  pulse  and  the 
S.I.T.  gate  was  varied  between  zero  and  one  hundred  jisec.  The  targets 
were  irradiated  at  a central  average  flux  level  of  5 x 10^  watts/cm2. 

Figure  10  shows  the  resulting  spectrum  between  3500  K and  5000  R 
with  a 1 p.sec  delay  between  laser  and  gate.  A black  body  spectrum,  with 
emissivity  and  temperature  as  free  parameters,  was  fitted  to  the  data  at 
7 points  using  the  least  mean  square  technique.  This  yielded  a temperature 
of  3707°K  _+  3%  and  an  emissivity  of  0.  70  +_  30%.  The  Chi  squared  test 
showed  a reasonable  fit.  The  data  was  also  fitted  in  the  least  mean  square 
sense  with  only  temperature  as  the  free  parameter  and  with  emissivity  held 
at  the  handbook  value  of  0.97.  This  gave  a temperature  of  3580°K  _+  3%. 
Along  with  the  experimental  spectrum,  Fig.  10  shows  a 3580°K  spectrum 
that  has  been  folded  with  the  measured  absolute  spectral  sensitivity  of  the 
spectrometer.  It  is  evident  that  the  observed  spectrum  can  be  attributed  to 
a hot  black  body  and  that  there  are  essentially  no  lines  due  to  superheated 
impurities. 
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Fig.  10  Carbon  Spectrum.  Dashed  line  is  a best  fit  3580°K  black  body 
spectrum.  Relative  timing  of  laser  and  spectrum  acquisition 
is  shown  in  lower  oscillogram 


Similar  results  were  obtained  for  a delay  of  4 fxscc  between  the 
start  of  the  laser  pulse  and  the  start  of  the  S.  I.  T.  gate  pulse.  The  tempera 
ture  had  increased  to  4066°K  _+  assuming  0.97  emissivity.  As  the  delay 
was  increased  further  to  7 jiSec,  the  emission  increased,  but  it  was  no 
longer  possible  to  make  a significant  fit  to  the  data  with  a black  body  curve. 
It  appeared  that  the  emissivity  towards  the  red  end  of  the  spectrum  was  low 
compared  to  the  emissivity  at  the  blue  end,  a feature  that  persisted  at 
longe  r delays 

Interferograms  and  shadowgrams  were  made  of  the  carbon  target  to 
show  its  vaporization  history  (Fig.  11).  A 300  nsec  Xe  pulsed  laser  at 
53  j3  R was  used  fc  r illumination.  The  pictures  were  made  on  a shot-to- 
shot  basis  on  fresh  target  spots,  but  reproducibility  was  good.  It  is  evident 
that,  lor  the  longer  delays,  the  spectrometer  is  viewing  the  target  surface 
through  a layer  of  blow-off  vapor.  If  the  vapor  and  the  target  have  differ- 
ent temperatures  and  the  optical  depth  of  the  vapor  is  less  than  unity  and 
varies  through  the  observed  spectral  regime,  non-black  body  profiles  would 
be  expected.  A useful  extension  of  the  present  data  would  include  observa- 
tions of  vapor  alone  by  sighting  the  spectrometer  at  right  angles  to  the 
target. 

Preliminary  emission  spectra  have  also  been  obtained  with  A1  tar- 
gets. Figure  12  shows  the  emission  in  the  wavelength  range  3500  to 
5000  R.  The  spectrometer  was  oriented  to  view  the  target  at  an  angle  of 
45°  as  with  the  carbon  targets.  For  Fig.  12  the  S.I  T.  gate  opened  coin- 
cidentally with  the  start  of  the  laser  pulse  and  it  remained  open  for  10  gisec. 
The  laser  flux  was  adjusted  to  a level  just  belov/  the  breakdown  threshold 
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Fig.  11  Carbon  Vapor  Blowoff.  Xe  laser  used  for  illumination. 
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Fig.  12  A1  Target  Emission  Spectrum  during  First  10  jisec  of  Las 
Irradiation 


Several  atomic  line  and  molecular  band  features  appear  superim- 
posed on  a broad  continuum.  In  particular,  the  neutral  A1  lines  at  3961  R 
and  3944  R are  quite  prominent.  The  molecular  band  system  is  due  to  AlO. 
At  this  point  it  is  not  clear  whether  the  AlO  is  due  to  oxidation  of  A1  vapor, 
or  decomposition  of  a surface  layer  of  hydroxide  AlO  (OH)  or  oxide  Al^. 

Some  shots  showed  the  as  yet  unidentified  line  at  « 4228  R as  the 
only  prominent  feature.  Candidates  are  singly  ionized  aluminum  at  4228  J? 
neutral  copper  at  4231  R and  neutral  nitrogen  at  4230  R.  Identification  is 
complicated  by  the  fact  that  wavelength  uncertainty  is  +_  3 J?  at  the  low  dis- 
persion employed.  Copper  is  a 4.  5%  by  weight  constituent  of  the  2024  alloy 
used  in  this  investigation.  The  boiling  point  of  copper,  2868°C.  is  higher 

than  Al,  2740°C,  so  it  would  be  surprising  if  this  material  left  the  surface 
preferentially. 

C.  DISCUSSION 

There  are  two  principal  questions  regarding  the  data  that  has  been 
presented.  The  first  question  deals  with  the  time  to  reach  vaporization 
temperature  at  a given  flux  and  the  second  asks  whether  a plasma  breakdown 
is  to  be  expected  in  the  observed  times  given  the  presence  of  vapor. 

The  time,  r,  required  for  the  surface  of  a uniform  half  space  of 

material  to  increase  in  temperature  by  AT  given  a heat  input  to  the  surface 
of  al  is 


* rp(c  at  +H)v 

■ 7 4— EJI — J 


mass  density,  Cp  is  the  specific 


where  y is  the  thermal  diffusivity,  p is  the 

W per  unit  ma.s,  H i.  the  heat  per  unit  mass  due  to  any  phase  change.  I 
ia  the  incident  power  per  unit  area  and  o ie  the  energy  absorptivity  of  the 


surface. 


Given  relatively  high  absorptivities  at  10.  6^,  relatively  low  dif- 
fusivities  and  low  vaporization  temperature  in  the  case  of  zinc,  the  theoret- 
ical vaporization  times  are  shorter  than  the  resolution  time  of  the  appa- 
ratus. Thus,  the  observation  of  prompt  vaporization  if.  not  inconsistent 
with  the  model  for  these  materials.  Similarly,  a carbon  surface  is  pre- 
dicted to  reach  the  observed  temperature  of  3580°K  in  time  of  1.1  x 106/ 
(Ql)  sec  using  a diffusivity  of  9.  3 x 10'^  cnZ/sec,  heat  capacity  of  0.  69 
joules/°K-gm  and  a density  of  1.7  3 gm/cm3.  The  flux  was  5 x 106  watts/ 
cm  and  the  absorptivity  is  expected  to  be  on  the  order  of  50%  leading  to  a 
predicted  time  of  less  than  1 psec.  Sublimation  of  the  surface,  which  can- 
not be  neglected  at  this  temperature,  would  tend  to  thermostat  the  surface 
temperature. 

In  the  case  of  aluminum  targets,  there  is  a large  discrepancy  be- 
tween the  observed  vaporization  time  of  3 psec  or  less  at  an  average  inci- 
dent flux  of  ~ 8 x 10  watts/cm2  and  the  theoretical  prediction  of  520  ^sec. 
This  latter  value  is  obtained  under  the  assumptions  of  3%  absorptivity,  a 
diffusivity  of  0.  35  cm  /sec,  and  an  energy  per  unit  volume  to  reach  vapor- 
ization, p(cpAT  + H),  of  10.  3 kj/cm3. 

A plausible  explanation  for  the  fast  vaporization  can  be  based  on  the 
assumption  of  surface  imperfections.  Consider  a surface  that  has  detached 
flakes  or  a rough  surface  such  as  shown  schematically  in  Fig.  13.  It  is  as- 
sumed that  the  scale  size  of  the  imperfections,  6,  is  less  than  the  thermal 
diffusion  distance,  V^t",  where  t is  the  time  for  the  surface  irregularity  to 
reach  vaporization  temperature.  By  equating  the  absorbed  energy  per  unit 
area  of  the  imperfection  divided  by  the  dimension  of  the  imperfection  to  the 
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energy  per  unit  volume  necessary  to  reach  vaporization  temperature,  a 
rough  criterion  for  the  vaporization  time  is  obtained 

p(c  AT  + H)  6 


t » 


al 


Alternately,  this  expression  may  be  solved  for  the  dimension  6 given  an 
observed  vaporizacion  time,  t. 

Using  the  observation  that  the  aluminum  targets  were  vaporizing  in 
3 fj  sec  or  less  at  8 x ' 06  watts/cm2,  the  inference  is  that  .mperfections  of 
dimens, on  less  than  on  micron  were  responsible.  This  dimension  is  less 
than  the  approximately  10  M diffusion  depth  at  3 |isec,  consistent  with  our 
original  assumption. 

The  evidence  for  surface  imperfections  is  multifold.  Time  inte- 
grated and  time  resolved  photographs  show  a number  of  localized  bright 
spots  on  the  surface  during  irradiation.  Scanning  electron  micrographs  of 
Al  surfaces,9  have  reve -led  surface  features  that  resemble  the  assumed 
flakes.  Moreover,  after  irradiation  with  a CO.,  laser  these  micrographs 
indicate  oxidation  of  the  flakes  as  evidenced  by  their  whiter  appearance 
which  would  result  from  enhanced  reflection  of  the  scanning  electrons  by  a 
nonconducting  feature.  Final'/,  there  is  the  observation  made  in  this  re- 
port that  vaporization  did  not  occur  when  the  surface  had  been  melted  by  a 
previous  laser  shot.  Such  melting  would,  of  course,  tend  to  remove  the 
micro -imperfections . 

The  time  to  heat  Al  vapor  to  full  ionization  as  a function  of  incident 
10.  6 |i  flux  has  been  calculated  by  Thomas  and  Musai,  and  by  Boni,  Cohen 
and  Su.6  Both  calculations  attribute  electron  heating  to  electron  neutral 
inverse  brems sf  rahlung  at  low  electron  density  and  electron  ion  inverse 
bremsstrahlung  at  high  electron  density.  Electrons  lose  energy  through 
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ionization  and  elastic  collisions  with  heavies.  There  are  differences 
between  the  calculations.  Thomas  and  Musal  assume  that  ionization  must 
take  place  between  the  ground  state  and  the  continuum,  while  Boni  et  al.  , 
argue  that  ionization  proceeds  through  transitions  to  the  first  resonance 
level.  For  constant  absorbed  laser  flux  these  assumptions  lead  to  different 
electron  temperatures  but  not  radically  different  breakdown  times.  The 
reason  is  that  absorbed  laser  power  is  approximately  balanced  by  power 
going  into  ionization  so  the  electron  temperature  will  rise  to  whatever  value 
is  necessary  to  give  the  requisite  ionization  rate.  (At  fluxes  much  below 

10  watts/cm2  the  calculations  also  differ  due  to  their  treatment  of  elastic 
cooling  losses . ) 

For  aluminum  vapor  leaving  the  target  surface  under  conditions  of 
temperature  and  pressure  corresponding  to  steady-state  vaporization  with 
an  incident  flux  of  107  watts /cm2,  Thomas  and  Musal  predict  breakdown 
(i*  G*  ’ fuU  ionization)  ^ 0.  1 jisec.  More  realistically  they  also  compute  a 
breakdown  with  initial  conditions  of  1 atmosphere  pressure  and  the  vapor  at 
the  Al  boiling  temperature  (~2750°K).  If  the  vapor  is  coming  from  small 
defect  points  on  the  surface  as  opposed  to  uniform  planar  vaporization,  the 
pressure  of  the  vapor  would  be  expected  to  rapidly  relax  to  1 atmosphere. 
The  predicted  breakdown  times  by  Thomas  et  al.  and  Boni  et  al.  , fall  in  the 
range  of  ten  to  fifty  microseconds  at  a flux  of  107  watts/cm2.  These  com- 
pare favorably  with  many  of  the  observed  breakdown  times.  However, 
some  breakdowns  occurred  in  shorter  times.  One  possible  explanation 
would  be  the  initiation  of  vaporization  at  more  sites  which  would  lead  to  a 
situation  corresponding  more  to  the  steady-state  initial  conditions  which 
are  calculated  to  have  shorter  breakdown  times. 
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III.  ENHANCED  TARGET  DAMAGE 


For  most  of  the  studies  involving  onset  of  target  vaporization,  the 
laser  flux  was  adjusted  to  be  on  the  borderline  of  the  threshold  for  break- 
down. Ihe  surfaces  showed  evidence  of  melt,  the  extent  of  which  increased 
with  laser  power.  Wherever  the  laser  power  was  increased  substantially 
above  the  breakdown  threshold,  on  the  order  of  10%  or  more,  gross 
cratering  of  the  surface  was  observed.  Figure  14  shows  the  change  in 
surface  appearance  of  2024  T4  A1  target  samples  as  the  laser  flux  is  in- 
creased. Each  photograph  shows  the  effect  of  a single  laser  pulse.  Break- 
down was  observed  on  all  but  the  lower  right  sample.  It  is  evident  that 
the  degree  of  cratering  increases  with  laser  flux. 

A.  DESCRIPTION  OF  EXPERIMENT 

The  samples  were  all  taken  from  the  same  sheet  of  Al  with  the 
surface  finish  in  an  "as  received"  condition.  Before  irradiation  the  sur- 
faces were  cleaned  with  acetone  followed  by  methyl  alcohol. 

Open  shutter  photographs  were  taken  to  show  the  extent  of  the 
breakdown  plume  and  a photodiode  was  trained  on  the  target  to  record 
the  time  history  of  light  from  the  breakdown. 

After  irradiation  the  samples  were  accurately  cross  sectioned  along 
the  line  forming  the  diameter  of  the  damage  spot.  Employing  standard 
metallur graphic  techniques,  the  samples  were  mounted  in  bakelite  and  the 
cross  sectioned  surfaces  were  polished  and  etched,  Fig.  15.  With  a 
measuring  microscope  the  depth  of  removed  material  on  axis,  and  the 
depth  ol  remaining  melted  material  on  axis  were  recorded,  Table  I. 
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Fig.  14  Cratering  on  A1  Targets  as  a Function  of  Incident  Energy. 

Surface  plasma  breakdown  occurred  on  all  but  the  18  joule  shot. 
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ENHANCED  COUPLING 
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Relative  instantaneous  laser  power  was  monitored  with  a gold  dooed 
germanium  photodiode  and  the  time  integrated  photodiode  signal  was 
standardized  against  total  energy  output  as  measured  with  a calorimeter. 
The  theoretical  value  for  the  fluence,  Frp,  (energy  per  unit  area)  on  axis 
at  the  focal  point  given  a circular  annular  aperture  with  output  coupling 
C is 
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T = 


(2.  44) 


2 n 
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E 
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where  DT  is  the  diameter  of  the  first  Airy  disk  at  the  focus  and  E is 
the  total  laser  output  energy.  DT  was  calculated  from  the  optical  system 
focal  length  and  measured  from  lucite  burn  patterns  as  0.  29  cm.  The 
output  coupling,  as  defined  by  the  geomefv  of  the  Breakdown  Laser's 
unstable  resonator  was  0.  80.  The  axial  fluences  employed  in  the  data 
reduction  to  be  presented  were  reduced  by  a factor  of  0.  72  from  ideal 
based  on  previous  (Ref.  11)  detailed  measurements  in  which  the  far-field 
intensity  distribution  was  mapped  by  moving  a narrow  slit  across  the 
focal  points. 

B.  ANALYSIS  OF  THERMAL  COUPLING 

An  effective  axial  thermal  coupling  coefficient  may  be  defined  by 
the  ratio 


pQ  (h  + h ) +pLh 
a_  m r m ^ s s 

F 

where  p is  the  density,  2.  7 gm/cm  , is  the  energy  per  unit  mass  to 

bring  the  material  up  to  and  including  the  molten  state  (1. 05  kj/gm) 
and  Lg  is  the  energy  to  bring  the  solid  up  to  but  not  including  the  melt 
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temperature  (651  j/gm).  The  depth  of  removed  material  is  h^,  the  depth 
of  remaining  solidified  melt  is  h , and  an  effective  thermal  penetration 
depth  is  given  as  hg.  The  total  energy  per  unit  area  incident  on  the  axis 
of  the  target  is  given  bv  F. 

The  thermal  penetration  distance  is  obtained  from  the  temperature 
distribution,  T(z),  in  the  solid  phase 


h 

s 


(T(z)  - Tg] 


where  T^  is  the  melt  temperature,  Tg  is  the  ambient  temperature  deep 
in  the  solid  and  z is  a coordinate  extending  from  the  melt- solid  boundary 
into  the  solid.  During  the  laser  pulse,  the  melt  boundary  is  penetrating 
into  the  solid  at  velocity  v which  is  assumed  to  be  a constant  equal  to  the 
final  penetration  distance  divided  by  the  laser  pulse  time,  T 


v = 


h + h 
r m 


An  equation  for  the  temperature  distribution  is  most  easily  written 
by  considering  the  thermal  diffusion  equation  in  a system  of  coordinates 
in  which  the  melt- solid-boundary  is  stationary. 


9T  9T  2 92T 

9t  ' v 9z  " X _ 2 
dz 

2 

The  diffusivity  of  the  solid  is  taken  as  a constant  (0.  5 cm  /sec  for  Al). 
This  equation  has  the  solution  satisfying  boundary  conditions 
T(t  = 0,  z)  = T , T(t,  z = 0)  = Tm  T(t,  z^  »)  = T , Ref.  12. 
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T-  T 
s_ 

T - T 
m s 


_1 

2 


, z+vt 
erfc + e 


- vz 

X 


erfc 


z-  vt 


2v/xt 


2^Xt 


where  erfc(x)  is  equal  to  one  minus  the  error  function,  Tg  is  the  temperature 
deep  in  the  solid,  and  Tm  is  the  melt  temperature.  The  space  integral 
between  z = 0 and  z = 00  may  be  done  with  the  help  of  relations  given  in 
Ref.  I 3. 

h = - erf  - ; + (Xt)1^  i1  erfc 

s v 2Vx7F  2V^ 

In  this  expression,  il  erfc(x)  is  the  first  integral  of  erfc(x)  between  the 
limits  of  x and  infinity.  This  function  is  tabulated  in  Ref.  13.  In  the  limit 
of  vanishing  velocity,  series  expansion  of  erf(x)  gives 

~-^<xt>,/2 


d v - 0 y/H 

the  well  known  result  for  thermal  diffusion  in  a stationary  medium.  In  the 
limit  of  velocity  greater  than  the  instantaneous  thermal  diffusion  velocity 


v » 


Vv/t 


a result  that  may  also  be  obtained  by  consideration  of  the  steady-state 
problem  with  6T/ 0t  = 0. 

The  computed  values  for  hg  with  t = Tp  are  listed  in  Table  I.  It  is 
noted  that  the  thermal  penetration  depths  are  all  significantly  less  than  both 
the  Airy  disk  diameter  and  the  diameter  of  the  melt  affected  zone.  This 
justifies  a one-dimensional  approach  to  the  diffusion  problem. 
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From  Table  I it  is  seen  that  the  effective  axial  coupling  increased 
from  3.  6%,  a value  that  is  not  unreasonable  for  CW  coupling,  to  7.  8%  in 
the  presence  of  a strong  breakdown.  Very  little  melt  remained  in  the  case 
of  the  strong  breakdown 
C.  MELT  REMOVAL 

Some  of  the  melt  was  undoubtedly  removed  in  the  form  of  droplets, 
since  trails  of  burning  droplets  were  observed  in  time  exposed  photographs 
and  splash  patterns  radiated  from  the  edge  of  the  crater.  Two  actions 
could  contribute  to  the  melt  removal.  A radial  pressure  gradient  is  created 
in  the  liquid  due  to  strong  vaporization  on  the  laser  axis  vs  little  or  no 
vaporization  as  the  radius  of  the  first  Airy  disk  is  approached.  The  radial 
pressure  gradient  would  accelerate  the  liquid  tangent  to  the  surface  in  a 
circularly  symmetric  pattern.  A simple  argument  may  be  made  to  show 
that  this  effect  is  only  to  be  expected  at  small  spot  sizes.  From  Bernouli' s 
equation  or  from  energy  arguments,  the  average  radial  fluid  velocity,  U, 
will  be  approximately  given  by 

pu2«  P 

7 

where  P is  the  on  axis  pressure.  For  A1  vaporizing  at  a flux  near  10 

2 7 2 

watts/cm  P may  be  shown  to  be  of  the  order  of  10  atm  (10  dyne/cm  ). 

In  order  for  the  fluid  to  accelerate  out  of  the  crater  within  a pulse  time, 

the  velocity  must  satisfy 

U > t/t 
P 

where  r is  the  laser  spot  size.  Thus 

r < Tp  VWp 
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For  a laser  pulse  of  100  Usee  this  implies  that  r must  be  less  than  a few 
millimeters.  Since  the  spot  diameter  employed  was  approximately  3 mm 
diameter  this  mechanism  is  not  inconsistent  with  observations. 

A second  mechanism,  which  is  not  limited  to  small  spots,  could  also 
explain  the  observed  cratering.  Associated  with  the  overpressure  at  the 
target  surface  due  to  vapor  blowoff  and  plasma  generation  there  is  a radial 
outflow  of  vapor  and  heated  air.  This  outflow  is  evident  in  the  interfero- 
grams  presented  in  Section  I.  Analysis  of  wave  generation  on  the  molten 
aluminum  "sea”  via  the  Kelvin-Helmholtz  instability  shows,  Ref.  14,  a 
half  wavelength,  X/2.  for  maximum  growth  of 


X _ ?_££  1 

2 = 2 p'  Uf2 


with  an  exponentiation  time  of 


11  oa 


where  0 is  the  molten  metal  surface  tension  (890  dyne/cm  for  All  and  U' 
is  the  velocity  of  the  vapor  of  density  p'  tangent  to  the  surface.  For 
velocities  near  sonic  in  air  (3.4  x 104  cm/sec,  and  ho.  high  pressure  gas 
(3000°K,  10  atm)  with  a density  of  1.  2 X 10‘3  gm/cm3,  the  half  wavelength 
and  growth  time  are  respectively  26p  and  2 psec.  Although  these  values 
indicate  strong  wave  formation  daring  the  laser  pulse  time  it  should  be 
noted  that  the  theory  totally  neglects  boundary  layer  effects  which  demand 
that  the  vapor  velocity  be  zero  at  the  surface  and  that  our  choice  of  vapor 
parameters,  if  not  unreasonable,  was  at  least  somewhat  cavalier.  Clearly, 
a more  complete  instability  theory  is  called  for  and  better  measurements 
or  calculations  of  vapor  flow  properties  above  the  surface  are  needed. 
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Once  waves  are  formed  on  the  surface,  a simple  argument  may  be 

made  as  to  whether  a droplet  spray  will  be  broken  away  from  the  wave 

peaks.  Consider  a protuberance  on  the  surface  with  characteristic 

dimensions,  D,  both  normal  and  parallel  to  the  surface.  The  aerodynamic 

2 2 

lift  forces  acting  on  the  protuberence  will  be  of  the  order  D p'U'  . When 
these  forces  are  greater  than  the  surface  tension  forces,  of  order  ctD, 
which  hold  the  protuberence  to  the  bulk  of  the  melt,  a free  droplet  will  be 
formed.  Thus,  projections  with  dimensions  satisfying 


jj  5" 

p'U*' 

will  form  a spray.  It  is  noted  that  this  critical  dimension  is  of  the  same 
order  as  the  half  wavelength  for  maximum  growth  of  the  Kelvin-Helmholtz 
instability.  Thus,  wave  generated  droplet  spray  remains  an  alternative 
for  explaining  the  observed  cratering.  This  mechanism  should  remain 
operative  for  large  target  spot  sizes  and  could  be  enhanced  by  air  flows 
on  the  order  of  Mach  1 over  the  target  surface. 
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